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Abstract: We propose an efficient approach to improve few-cycle soHton 
compression with cascaded quadratic nonfinearities by using an engineered 
multi- section structure of the nonUnear crystal. By exploiting engineering of 
the cascaded quadratic nonlinearities, in each section soliton compression 
with a low effective order is realized, and high-quality few-cycle pulses with 
large compression factors are feasible. Each subsequent section is designed 
so that the compressed pulse exiting the previous section experiences 
an overall effective self-defocusing cubic nonlinearity corresponding to 
a modest soliton order, which is kept larger than unity to ensure further 
compression. This is done by increasing the cascaded quadratic nonlinearity 
in the new section with an engineered reduced residual phase mismatch. 
The low soliton orders in each section ensure excellent pulse quality and 
high efficiency. Numerical results show that compressed pulses with less 
than three-cycle duration can be achieved even when the compression 
factor is very large, and in contrast to standard soliton compression, these 
compressed pulses have minimal pedestal and high quality factor. 

© 2012 Optical Society of America 

OCIS codes: (190.5530) Pulse propagation and temporal solitons; (190.7110) Ultrafast non- 
linear optics; (320.5520) Pulse compression. 



References and links 

1. E. Seise, A. Klenke, S. Breitkopf, M. Plotner, J. Limpert, and A. Tiinnermann, "Coherently combined fiber laser 
system delivering 120 J femtosecond pulses," Opt. Lett. 36, 439-441 (2011). 

2. M. Nisoli, S. De Silvestri, and O. Svelto, "Generation of high energy 10 fs pulses by a new pulse compression 
technique," Appl. Phys. Lett. 68, 2793-2795 (1996). 

3. CP. Hauri, W. Kornelis, F. W. Helbing, A. Heinrich, A. Couairon, A. Mysyrowicz, J. Biegert, and U. Keller, 
"Generation of intense, carrier-envelope phase-locked few-cycle laser pulses through filamentation," Appl. Phys. 
B 79, 673-677 (2004). 

4. X. Liu, L. Qian, and F. W. Wise, "High-energy pulse compression by use of negative phase shifts produced by 
the cascaded x^'^^ ■ X^^^ nonlinearity," Opt. Lett. 24, \111-\119 (1999). 

5. S. Ashihara, J. Nishina, T. Shimura, and K. Kuroda, "Soliton compression of femtosecond pulses in quadratic 
media," J. Opt. Soc. Am. B 19, 2505-2510 (2002). 

6. S. Ashihara, T. Shimura, K. Kuroda, N. E. Yu, S. Kurimura, K. Kitamura, M. Cha, and T. Taira, "Optical pulse 
compression using cascaded quadratic nonlinearities in periodically poled lithium niobate," Appl. Phys. Lett. 84, 
1055-1057 (2004). 

7. X. Zeng, S. Ashihara, X. Chen, T. Shimura, and K. Kuroda, "Two-color pulse compression in aperiodically poled 
lithium niobate," Opt. Commun. 281, 4499-4503 (2008). 



8. J. Moses and F. W. Wise, "Soliton compression in quadratic media: high-energy few-cycle pulses with a 
frequency-doubling crystal," Opt. Lett. 31, 1881-1883 (2006). 

9. G. Xie, D. Zhang, L. Qian, H. Zhu, and D. Tang, "Multi-stage pulse compression by use of cascaded quadratic 
nonlinearity," Opt. Commun. 273, 207-213 (2007). 

10. X. Zeng, S. Ashihara, N. Fujioka, T. Shimura, and K. Kuroda, "Adiabatic compression of quadratic temporal 
solitons in aperiodic quasi-phase-matching gratings," Opt. Express 14, 9358-9370 (2006). 

11. X. Zeng, S. Ashihara, Z. Wang, T. Wang, Y. Chen, and M. Cha, "Excitation of two-colored temporal solitons in 
a segmented quasi-phase-matching structure," Opt. Express 17, 16877-16884 (2009). 

12. M. Bache, O. Bang, J. Moses, and F. W. Wise, "Nonlocal explanation of stationary and nonstationary regimes in 
cascaded soliton pulse compression," Opt. Lett. 32, 2490-2492 (2007). 

13. B. B. Zhou, A. Chong, F. W. Wise and M. Bache, "Few-cycle solitons in short strongly phase-mismatched 
frequency conversion crystals," Phys. Rev. Lett. 109, 043902 (2012), arXiv: 1109.4261. 

14. M. Bache, J. Moses, and F. W. Wise, "Scaling laws for soliton pulse compression by cascaded quadratic nonlin- 
earities," J. Opt. Soc. Am. B 24, 2752-2762 (2007). 

15. M. Bache, O. Bang, W. Krolikowski, J. Moses, and F. W. Wise, "Limits to compression with cascaded quadratic 
soliton compressors," Opt. Express 16, 3273-3287 (2008). 

16. J. Moses and F. W. Wise, "Controllable self-steepening of ultrashort pulses in quadratic nonlinear media," Phys. 
Rev. Lett. 97, 073903 (2006). See also arXiv:physics/0604170. 

17. R. Schick and T. Pertsch, "Absolute measurement of the quadratic nonlinear susceptibility of lithium niobate in 
waveguides," Opt. Mater. Express 2, 126-139 (2012). 

18. O. Gayer, Z. Sacks, E. Galun, and A. Arie, "Temperature and wavelength dependent refractive index equations 
for MgO-doped congruent and stoichiometric LiNbOs," Appl. Phys. B 91, 343 (2008). 

19. M. Bache, O. Bang, B. B. Zhou, J. Moses and F. W. Wise, "Optical Cherenkov radiation in ultra-fast cascaded 
second-harmonic generation," Phys. Rev. A 82, 063806 (2010). 

20. C. R. Phillips, C. Langrock, J. S. Pelc, M. M. Fejer, L Hartl, and M. E. Fermann, "Supercontinuum generation in 
quasi-phase-matched waveguides," Opt. Express 19, 18754-18773 (2011). 

21. M. Bache, O. Bang, B. B. Zhou, J. Moses, and F. W. Wise, "Optical Cherenkov radiation by cascaded nonlinear 
interaction: an efficient source of energetic few-cycle near- to mid-IR pulses," Opt. Express 19, 22557-22562 
(2011). 

22. M. Bache, R. Schick, "Review of measurements of Kerr nonlinearities in lithium niobate: the role of the delayed 
Raman response," arXiv: 121 1.1721 (2012), http://arxiv.org/abs/1211.1721. 

23. H. Guo, X. Zeng, B. B. Zhou, and M. Bache, "Electric-field modeling and self-steepening counterbalance of 
cascading nonlinear soliton pulse compression," arXiv: 1210.5903 (2012), http://arxiv.org/abs/1210.5903. 

24. M. Bache, H. Guo, B. B. Zhou and X. Zeng, "The anisotropic Kerr nonlinear refractive index of the beta-barium 
borate (j8-BaB204) nonlinear crystal," arXiv: 1209.3 158 (2012), http://arxiv.org/abs/1209.3158. 

25. C. Myatt, N. Traggis, and K. L. Dessau, "Optical fabrication: Optical contacting grows more robust," Laser Focus 
World 42, 95-98 (2006). 



1. Introduction 

High-energy pulsed lasers by employing chirped-pulse-amplification (CPA) technique are de- 
veloped rapidly both in bulk solid-state and fiber amplifiers. Recently a power scaling approach 
by using coherent combination of several channels of femtosecond fiber CPA systems was 
demonstrated to hugely increase the deliverable power [1]. Besides the popular Ti: sapphire 
technology, currently there is a lot of focus on CPAs at the wavelength of 1.03 fim by us- 
ing Yb-doped or 1.56 /im by using Er-doped gain materials. However, their pulse durations 
are severely limited to several hundreds of femtosecond due to the narrow gain bandwidth of 
the active materials. Generation of optical few-cycle pulses inside the CPA systems requires 
nonlinear phase shift to broaden the spectrum, which may lead to the distorted pulses due to 
excessive detrimental nonlinear effects. 

Pulse compression to few-cycle duration from hundreds of femtosecond duration out of the 
energetic solid-state and fiber amplifiers, would give a stable and compact optical source for e.g. 
ultrafast pump-probe spectroscopy and biology. One way of achieving this would be to use gas- 
filled cells [^] or hollow capillaries [ ] to achieve spectral broadening followed by dispersive 
optics like prisms or gratings to compress the pulses temporally, but these techniques, while 
standard, suffer from being hard to control and being very bulky. 

An alternative is to use soliton compression from self-defocusing cascaded nonlinearities 



(e.g., phase-mismatched second-harmonic generation, SHG): when they are strong enough to 
overcome the inherent cubic (Kerr) self-focusing nonUnearity in the nonUnear materials, an ef- 
fective self-defocusing nonlinearity is formed [4] without beam filamentation or break-up of 
the transverse part of the beam, problems often encountered with self-focusing nonlinearities, 
and the compressible pulse energy is in principle only limited by the aperture of the nonlinear 
medium. Furthermore solitons can be generated in the visible and near-IR regimes [^ ]. This is 
possible because the effective negative nonlinear refractive index can be counterbalanced by the 
positive (normal) dispersion to create a soliton. When taking the effective soliton order from 
the total defocusing nonlinearity substantially larger than unity, in the initial propagation stage 
the input pulse will experience a dramatic self-compression effect. This soliton compression 
works in the same way as standard soliton compression with Kerr nonlinearities in fibers. Soli- 
ton compression by using cascaded quadratic nonlinearities was experimentally observed from 
120 fs down to around 30 fs at the wavelength of 800 nm [ ] and at telecom wavelengths [ , ]. 
12 fs (3 optical cycles) compressed from 1 10 fs was achieved at 1250 nm [ ]. Multi-stage pulse 
compression in sequence by use of cascaded quadratic nonlinearity was proposed to achieve 
extreme compression factor from several ps to 30 fs [ ]. Quasi-phase-matching (QPM) tech- 
nique offer additional degrees of freedom controlling cascaded nonlinearities, depending on the 
effective phase mismatch. Through engineering the QPM structure, adiabatic and high-quality 
pulse formation of quadratic soliton was presented [10, 1 1]. 

Recently a temporal nonlocal model gave insight in optimal phase-mismatch ranges, where 
few-cycle pulse durations are possible [ ] and this led to an experimental demonstration of 17 
fs compressed soliton (less than 4 optical cycles) from 47 fs near-IR pulses through cascaded 
quadratic soliton compression (CQSC) in a short standard lithium niobate crystal [ ]. 
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Fig. 1. A schematic of multi- section structure of nonlinear crystals, (a) Cascaded nonlinear- 
ity increases and (b) the effective phase mismatch deceases upon propagation by controlling 
the local domain period Aj in (c) QPM grating or the angle 6j between the FW input and 
the optical z axis of the crystal in (d) bonded BBO crystals. 



However, the soliton compression mechanism used in the connection with cascaded non- 
linearities has similar drawbacks to the higher-order soliton compression known from optical 
fibers: the larger the soliton number, the larger the unwanted pedestal, and the less energy re- 
mains in the compressed spike. Such a feature is detrimental to practical applications, and is 



in particular a problem with large compression factors, and this is exactly the challenge faced 
here where we seek to compress longer >> 100 fs pulses to sub-20 fs duration. 

Here we propose an efficient approach to improve few-cycle soliton compression by en- 
gineering multi-section structures in nonlinear crystals, such as quasi-phase-matching (QPM) 
structures with reduced residential phase mismatch or bonded j3 -barium borate (BBO) crystals 
with different orientation of optical axes, as shown in Fig. 1. The effective phase mismatch 
(A^eff) in consecutive sections of nonlinear crystals decreases and always keeps A^eff > to 
make the cascading nonlinearity self-defocusing, and in this way the cascaded quadratic non- 
linearity is engineered through the phase mismatch, a property which is not possible in standard 
Kerr media. The basic idea is that each section is designed to make the effective soliton order 
A^eff > 1 by adjusting the residual (or QPM) phase mismatch, and so that at the exit of each sec- 
tion soliton compression is achieved. Such a compressed soliton will at this stage be in balance, 
so no further compression occurs unless we change the conditions, which is exactly what we 
propose. When this compressed pulse enters the new section the nonlinearity is engineered to 
make further compression occur. At the same time, each section is designed to keep a low effec- 
tive soliton number, as soliton compression of pulses with small soliton numbers can achieve a 
small pedestal and high pulse quality. The numerical results turn out to work remarkably well 
already as a two- section and three- section design. 

2. Propagation equations and nonlinear Schrodinger equation in the cascading limit 

The temporal dynamics of ultra-short pulses propagating in the nonlinear media, specifically 
BBO and LN, can be described by the coupled wave equations under the slowly varying enve- 
lope approximation [14, 15]: 



= (1) 
= (2) 



where Ej{z^t) is the complex amplitude of the electric field, pj and Gj involve quadratic and 
cubic nonlinearities, respectively, as pj{z) = (Oid^f^iz) / cnj and Oj = 3()L)jRe(;^(^))/8c^j. Jeff(^) 
is the quadratic nonlinear susceptibility along the propagation. The subscripts 1 and 2 corre- 
spond to the fundamental wave (FW) and second harmonic (SH) pulses. Here the phase mis- 
match Ako = k2 — 2k\ and the material dispersion coefficients k^^^ = d^kj/ dco^ \ (o=a)j ^ where 
kj = rijCOj/ c is the wave numbers and Uj is the refractive indices. Time t is measured in a frame 
of reference travelling with the group velocity of the FW pulse. The group velocity mismatch 
(GVM) is di2 = k^l^ - = V-} - v'l and Dj = 1^=2 '^.^f^TP^ is total dispersion opera- 
tor, which in frequency domain can be described as Dj{0) = kj{cOj -\- Q,) — Q,k[^'^ — kj{cOj). 
The cross-phase modulation (XPM) coefficient is T] = 2 (or 2/3) for type-0 (or type-I) SHG 
geometry. 

Under the conditions of the cascading limit (AkL ^1) and the negligible cross-phase- 
modulation (XPM) contribution, the soliton compression based on cascaded phase mismatched 
SHG process can be easily understood by reducing Eqs. (1-2) to a normalized nonlinear 
Schrodinger equation (NLSE) for the FW field Ui (see details in Ref. [15]): 

{i^ + ^ - (nLc - NlM P = iNino^m (3) 
The dimensionless field Ui is normalized by the peak of the input electric field Ei^ and 



the time. The propagating distance are rescaled di^ T = {t — z) /Tq and = z/Ld, where 
= Tq /\k^p\ and Tq are the FW dispersion length and input pulse duration. The dimen- 

sionless phase mismatch and GVM is Aj8 = A^q^j and 5 = Ji27o/|^P^|. Here we assume 
that A/:eff > so the cascading nonlinearity is self-defocusing, and that the FW GVD is pos- 
itive (normal dispersion). The effective soliton number is defined as N^^ = N^^^^ — N^^^ = 
Lj()L)i/o( Incase I ~^Kerr)/^' where Iq is input intensity and A/casc and A^Kerr are the dimensionless 
quadratic and cubic soliton numbers and ^^asc ^ ~<^eff/^^eff is the Kerr-like cascading and n^^^^^ 
is the self-focusing material Kerr nonlinearities, respectively. 

Equation (3) shows that with cascaded SHG, solitons can be excited in much the same way 
as in a Kerr medium [ ] : In order for few-cycle pulse compression to occur based on phase- 
mismatched SHG, an effective self-defocusing nonlinearity A/eff > 1 with normal dispersion is 
required. In the same way as for solitons in a pure Kerr medium, pulse compression with higher- 
order solitons suffer from low pulse quality and a large unwanted uncompressed pedestal. The 
right-hand side of Eq. (3) contains the first-order contribution (in the weakly nonlocal approx- 
imation) of the cascading response. It is a GVM-induced contribution, which acts in a similar 
way as self- steepening [16]. 

3. Numerical results and discussions 

To obtain accurate numerical simulation of few-cycle pulse region, we use the full coupled SHG 
propagation equations under the slowly-evolving-wave approximation including exact disper- 
sion, self- steepening and electronic Kerr effects [ ] (see [ ] for more details). We consider 
5% MgO doped periodically poled lithium niobate (PPLN: MgO) for pulse compression at the 
wavelength of 1.56 /im, where high-energy Er-doped fiber amplifiers operate. Type-0 phase- 
mismatched SHG is used for soliton compression because of its large quadratic nonlinear com- 
ponent (at 1520 nm the value J33 = 20.6 pmA^ was measured recently [ ] and the electronic 
Kerr nonlinear refractive index is chosen to be ^^Kerr = 3Re{x^^^)/4nleoc = 23 • lO-^^mYw 
[13]). The refractive indices and the material dispersion coefficients are derived from Sell- 
meier's equation for MgO: LN [18]. The balancing point between quadratic and cubic non- 
linearities where |^casc| — ^Kerr determines the upper limit of the effective phase mismatch 
for soliton formation in normal dispersion, which is calculated to be less than 212 mm~^ in 
QPM and the effective phase mismatch is A^j,qpm = A^o — 2;r/ and A^ is the QPM period. 
Therefore we always choose A^i qpm in the first QPM section to be 180 mm~^ in the following 
simulations of soliton compression in multi-section QPM gratings. 

The input transform-limited FF pulse duration is 250 fs at full width at half-maximum 
(FWHM) at the pump wavelength of 1.56 /im. Figure 2(a)-(b) show the numerical simulations 
of soliton compression with different effective soliton numbers (A^i,eff) in single- section QPM 
structures. The pulse quality of the compressed FW is well maintained by keeping a low soliton 
order but the drawback is a limited smaller compression factor. Figure 2(c)-(d) shows the FW 
temporal and spectral evolution in a two-section QPM structure. The effective soliton number 
in the first QPM section A^i^eff is chosen to 3, corresponding to a peak intensity of 25 GW/cm^ 
when A^i QPM = 180 mm~^ The temporal propagation dynamics for such a configuration is 
shown in Fig. 2(b). The compressed pulse reaches its minimal duration of 34 fs after propa- 
gating the distance of Li = 0.4 Li j (the dispersion length is Li j = 211 mm in the first QPM 
section) and continues to compress to 14 fs in the second QPM section, which has A/:2,qpm = 
140 mm~\ which corresponds to A^2,eff =1.6 for a 34 fs pulse. In the spectrum, strong broad- 
ening of the FW is observed, and the soliton formation also leads to a soliton-induced optical 
Cherenkov wave around the wavelength of 3 /im, as shown in Fig. 2(d). The FW pulse has 
a slight oscillation upon propagation in the onset of the second QPM and then continues to 




Fig. 2. Numerical simulation of compressed pulses of Tin = 250 fs FWHM at the pump 
wavelength of 1.56 jim in one- and two-section QPM structures. FW temporal evolutions 
of (a) A^i^eff =10^ ^^QPM = 150 mm~^ and (b) A^i^eff = 3, A^qpm =180 mm~^ in one- 
section QPM structure; (c) and (d) FW temporal and spectral evolution in two- section QPM 
structure with A/^i gff = 3, A^i^qpm = 180 mm~^ (A^2,eff =1-6, A^2,qpm = 140 mm~^) in the 
first (second) QPM. (e) and (f) Normalized intensities of the input, the optimal compression 
positions [white dashed lines, cuts in (a), (b) and (c)]. (g) Pulse durations and the ratio of 
peak intensities along propagation in two- section QPM. 

compress further and the peak intensity increases quickly, as shown in Fig. 2(f). 

Note that the dispersion length in the jth-section QPM is defined as Lj^d = ^fm/ l^f^l 
in is the pulse duration at the input of jth-section QPM. The second distance L2 for optimal 
compression is much shorter than Li , but it corresponds to 1 .4 L2,j and the dispersion length 
L2^d is 4 mm. The quality factor (Qc) of the compressed pulse, defined as the ratio of the 
fractional amount of energy carried by the central spike of the FW pulse over the input pulse 
energy, is as high as 0.63. The pulse compression of higher-order soliton compression (A/gff = 
10) in a single- section QPM is also shown in Fig. 2(e) for comparison (Qc is 0.22). The pulse 
quality is much worse than with two-section QPM: a large pedestal and trailing oscillations 
are observed, even though the main spike has pulse duration of 14 fs. The length required for 
multi- section QPM is in the example shown quite long (around 85 mm), but it can be reduced 
by increasing the soliton number. For instance, the total length required for pulse compression 
of the input FF pulse (250 fs) with the soliton number A/^i eff = 5 is 47 mm. This comes at a price, 
though, as the quality factor is reduced to 0.45 and the input intensity increases. Note also if 
LN is used at shorter wavelengths, the dispersion length becomes shorter due to the increase 
of material GVD. On the other hand, input Gaussian beam with a longer focus length than the 
device length is required and leads to the reduced pump intensity. Fortunately, our compression 
scheme is accomplished under a low soliton order, in which the pump intensity is low. Hence a 
laser pulse with an energy on the order of ten micro-Joules can meet the requirement. 

Figure 3(a) shows the simulations of soliton compression by changing the effective wave 
mismatch A^2,qpm of the second QPM section. The optimal pulse durations and peak intensities 
for different input pulse durations and soliton numbers are shown in Fig. 3(b). It is clearly seen 
that all of the best compression cases (i.e. pulses with the highest intensity and shortest duration) 
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Fig. 3. (a) Optimal compressed pulse intensities for selected values of the effective phase 
mismatch A^2,QPM of the second QPM and (b) pulse durations and peak intensities versus 
A^2,QPM- (c) Wavelength dependence of Ako from material dispersion, the resonant thresh- 

(2) 

old Ak^Y , the upper limit of A^c,LN for soliton in unpoled LN. (d) Compression windows 
versus the wavelength between the upper limit of A/:c,qpm (black line), the resonant thresh- 
old A^sr (red line) including full SH dispersion compression and the breakeven QPM value 
A^even,QPM (purple line) to achieve the same fi^asc unpoled LN. 



occur in second- section QPM when A^2,qpm approaches to 140mm~\ at which the critical 
transition occurs from the nonstationary to stationary regions in the nonlocal model [12], i.e., 
the nonlocal response is localized in the stationary regime and oscillating in the nonstationary 
region. 

The detailed discussions on the transition to the nonstationary regimes of BBO crystal were 
presented based on the nonlocal theory [12, 14, 15] . Here we show the nonlocal analysis of 
stationary and nonstationary regimes in type-0 phase-mismatching SHG geometry of PPLN: 
MgO in the normal dispersion regime. Now including the effective phase mismatch in QPM, 
the nonlocal response function in frequency domain is derived as 

R{^)=AkQPM/Mni{^) (4) 

and 

Aknl{^)=D2{^)-di2^^^kQPM=k2{C02^^)-2k^{C0i) - - k[^^ ^ (5) 

It is seen that when A^qpm is below a certain threshold A^qpm < A^sr, becomes reso- 
nant spectrally and consequently the SH spectrum develops a characteristic sharp resonance 
peak (in the time domain the SH pulse becomes stretched and can eventually span many pi- 
coseconds) [ ]. This is of course an unwanted property for efficient femtosecond pulse in- 
teraction, where a large bandwidth is required; when A^qpm is significantly larger than the 
threshold A^qpm > A/^sr, R{^) is non-resonant and holds octave- spanning bandwidths [ ]. 
Therefore A^qpm > A/:sr indicates the boundary to the stationary regime, where ultra-broadband 



non-resonant interaction is supported. If only including up to second-order SH dispersion, we 
can explicitly obtain the threshold A^s^^ = ^^2/2^2^^ t)y using full SH dispersion, A^^j. ^ 
143 mm~\ which is quite close to the numerical results of the phase-mismatch value of the 
optimal compression in the second QPM gratings. 

QPM reduces the effective quadratic nonlinearity Jgff (with at least 2/71 ) and the breakeven 
phase-mismatch value A^even,QPM = 4A^o /^^ ^ 130 mm~^ to achieve the same nonlinearity as 
in unpoled LN (^casc(^^o) = ^Lsc(^^even,QPM) = 38 • 10~^^m^/W). This phase mismatch value 
turns out to be in the non- stationary regime, while the inherent (material) phase mismatch Ako 
is large enough to be in the stationary regime and still be below the upper limit of A/:c,ln 
for soliton compression, as shown in Fig. 3(c). The cascading nonlinearity using QPM when 
A/:qpm > A/:sr, is therefore smaller than that in unpoled LN, and this property exhibits the 
advantage by using noncritical SHG for soliton compression in unpoled LN over QPM LN 
[13]. However in the multi-section structure, unpoled LN cannot be used for the first- stage 
compression. As in the following QPM section ^^^g^(A/:2,QPM) need to be increased, there is 
no QPM phase-mismatch space for second-stage compression because in the second stage, the 
cascaded nonlinear index needs to be larger than that in the first stage, seen in Fig. 1. One 
reason is that the breakeven phase mismatch value A^even,QPM in unpoled LN is close to the 
non- stationary regime and another is that the larger cascaded nonlinearity of second section 
PPLN requires smaller phase mismatching value because QPM modulation gives a pre-factor 
(2/71 from the first order). 

QPM can engineer In^asd <^eff/^^QPM by changing the residual phase mismatch with dif- 
ferent QPM domain lengths, even if this means sacrificing the cascading strength. Taking into 
account that optimal soliton compression occurs in the stationary regime (A^gff > A^sr) and that 
there is an upper limit of the effective phase mismatch for soliton formation in normal disper- 
sion, we obtain the compression window of A/:qpm in type-0 phase-matching SHG geometry 
of PPLN: MgO, as shown in Fig. 3(d). This shows enough A/:qpm space for fully exploit- 
ing quadratic nonlinearities by engineering the QPM structure, for example, ^^^^^(A^qpm = 
180mm-i) = 27 • IQ-^^mVw and ^^..(Ay^QPM = 140mm-i) = 35 • IQ-^^m^W. It is noted 
that the critical wavelength where the compression window opens is around 1.37 /im due to 
the reduced cascading nonlinearity, which longer than that in unpoled LN, where the window 
opens already at 0.9 /im. 

Figure 4(a)-(b) compare in one- and multi-section QPM design the quality factor and the 
optimum pulse duration versus the effective soliton number. In single- section QPM, a high 
soliton order is required to compress the pulse towards few-cycle duration, but this leads to 
much lower Qc; the compressed pulses with duration of an empirical limit of 13 fs and good 
quality factors (Qc > 0.4) are achieved in two- section QPM by using the input pulses of small 
soliton numbers (A/gff = 3, 4 and 5). The minimal pulse durations of 18 fs with a high quality 
factor 0.8 are achieved in three- and four-section QPM, and evidently the trend is that many 
sections and low soliton orders are favorable. The caveat is that the total propagation length 
becomes very long, which makes it difficult to be realized practically in a single QPM crystal, 
so as always one must find a good compromise in a final design. 

In fact, input pulses with a large soliton number will experience a stronger competing Kerr 
XPM effect, which is also self-focusing and reduces the overall self-defocusing nonlinearity. 
The effective soliton number corrected by XPM effect can be expressed as 

2 22 2 22 

■^efr(XPM) — ^casc ~ ^Kerr(SPM) ~ ^Kerr(XPM) — ^efl ~ ■^Kerr(XPM) 
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Fig. 4. Results of numerical simulations showing the optimum compression parameters 
versus A^gff i^i multi- section QPM gratings for (a) the minimum pulse duration, (b) quality 
factor. The filled mark with round symbols are the compression results in one section QPM 
with good quality factor (Qc > 0.4) and the square, star and diamond symbols represents the 
results in multi- section QPM (two, three and four) with high quality factor while durations 
are close to the empirical limit of 13 fs (green dash line). In (a), the arrow shows stronger 
pulse compression with the decrease of N^^^ and the black dash line is fit to the minimum 
duration available in one- section QPM gratings. The green dashed line in (a) denotes an 
empirical limit of 13 fs that was observed in the simulations. 



Thus the critical value of the effective soliton number for soliton behavior requires 

N,s>Ks=^Jl+^J^l (8) 



n^c^ j (^^Kerr Incase 1^1 ^i^f) • ^in IS the energy fluence of the input 

pulse, <l>in = 27in/in for a sech^-shaped pulse. A^^^^ is clearly dependent on the input fluence 
and the phase mismatch. In Fig. 4(a), we show that by using large soliton number stronger 
pulse compression occurs as A/^^^ decreases. 

One may question whether the approach works if the input parameters are different from the 
pre-designed values. The numerical results confirm the capability of generating high-quality 
compressed pulses, as shown in Fig. 5(a)-(c), when soliton compression in the first-stage QPM 
occurs early or later, which is beneficial from low-order soliton propagation by controlling 
quadratic nonlinearities. Figure 5(d) plots the peak intensifies of FW along the propagation 
direction as function of the first QPM length L\ . Soliton compression to few-cycle pulses still 
occur in second- section QPM even when L\ has become shorter or longer than the optimal 
length. Similarly, the capability to survive some fluctuation of input pulses is also seen in Fig. 
5(e). The final output pulses have slight degradation when 10% fluctuations of pulse intensity 
are inputted in the same QPM structure (L\ =44.5 mm and L2 = 3 mm). 

Recently Raman nonlinear contribution of lithium niobate in cascaded nonlinear interaction 
have been carefully considered in supercontinuum generation [20], few-cycle soliton compres- 
sion and optical Cherenkov radiation [13, 19, 21]. In the time scale of few-cycle pulses, Raman 
Kerr response of lithium niobate becomes quite complicated and relies on the property of ma- 
terial composition [22]. In the present study, Raman fractions = 0.14 (0.50) are used in the 
simulations, while keeping the same electronic Kerr contributions. 

In the case of small Raman fraction = 0.14, Fig. 6(a) shows a simulation in a single-section 
QPM structure. Compared with Fig. 2(b) where Raman was not included, the simulation shows 
that the pulse is sped up instead of being delayed. This is because in the normal-dispersion 



where <l>^ = 4/4 ^^^A^ 
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Fig. 5. Temporal evolutions of FW pulse (250 fs FWHM) in two-section QPM structures 
(^i,eff = 5): the length of first-section QPM is (a) 39.2 mm (b) 44.5 mm (optimal compres- 
sion position) and (c) 48.5 mm. (d) Peak intensifies of FW along propagation as function 
of the length of first QPM. (e) Comparison of the final output pulses with 10% fluctuation 
of input intensities in same QPM structure (Li = 44.5 mm and L2 = 3 mm). 
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Fig. 6. Numerical results of temporal evolution of FW pulse (250 fs FWHM) including 
Raman contribution, (a) one-section (b) two-section QPM when /r = 0.14 and (c) one- 
section QPM when /r = 0.50. 



regime, the soliton that is becoming red-shifted spectrally by the Raman effect will travel faster 
and it even overcomes the initial pulse delay from the GVM-induced self- steepening (right- 
hand side of Eq. (3)). For this value of the Raman fraction we can still implement the soliton 
compression to achieve high-quality few-cycle pulses by using multi-section QPM structure, as 
shown in Fig. 6(b). As the design imposes a decrease of A^qpm in the second- section QPM, the 
GVM-induced self- steepening increases (as it is proportional to 1/A/:qpm) to reduce the overall 
Raman influence and even leads to possible balance; this will discussed further in a separate 
publication [23]. However, when the Raman fraction is increased to fR = 0.50, a sort of modu- 
lational instability occurs as seen in Fig. 6(c). Note that it seems not to be soliton fission, which 
usually happens after the first compression point, as it seems that the pulse does not experience 
complete compression before pulse splitting happens. Therefore the Raman response of the 
nonlinear materials plays important role in the cascaded soliton compression and the dynamics 
will be published elsewhere. 
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Fig. 7. Evolution of the temporal profiles of FW pulses (Tin = 200 fs FWHM) in BBO 
crystals. Propagation in one section (a) A^eff = 3.5, A/^gff = 75 mm~^ (b) A/gff = 8, A^eff 
= 60 mm~^; (c) A^i = 75 mm~^ (A^i,eff = 3.5) and Ak2 = 58 mm~^ in two-section BBO 
crystal (Li = 93.9 mm and L2 = 4.5 mm); (d) pulse intensities at the optimal compression 
positions (dashed lines) in one-section (FWi) and two-section (FW2,i: first compression; 
FW2,2- second compression) BBO crystals. 

Finally we consider BBO as the nonlinear medium in type-I birefringent phase-matching 
SHG geometry, in which Raman contribution is negligible under the relevant experimental 
conditions. BBO has a decent quadratic nonlinear strength ( Jgff = 2.06 pmA^) and a very 
low cubic nonlinear refractive index ( ^^Kerr = 5.90- 10~^^m^/W) at the wavelength of 0.8 
fim [24]. Compared to LN, the FW GVD in BBO crystal is smaller, which leads to longer 
dispersion lengths. Therefore we use in the simulations shorter durations (200 fs FWHM) of 
input pulse and the wavelength of 1.030 /im, which is current state of the art for Yb-based 
solid-state amplifiers. As shown in Fig. 1(d), the phase mismatch can be altered by changing 
the angle 6 between the FW input and the optical z axis of the crystal. A multi-section structure 
of BBO crystals with different angles 6 may be constructed through some epoxy-free optical 
bonding techniques, as seen in [ ]. Alternatively, two BBO crystals could be put in sequence 
(one long and one short) where individual angle-tuning of the phase-mismatch values gives a 
large operational flexibility to achieve optimal pulse compression. Achieving huge compression 
factors from initial 30-ps to 30 fs pulses was proposed based on multi-stage pulse compression 
by use of cascaded quadratic nonlinearity in BBO [9], where additional beam expanders are 
needed to reduce the peak intensities. The difference to the current work is that we always keep 
soliton compression in the stationary regimes to achieve octave- spanning bandwidths that are 



large enough for supporting few-cycle pulses and use small soliton numbers in each section 
compression to produce high-quality pulses. 

The numerical results shown in Fig. 7 confirm that pulse compression in two-section BBO 
crystals is greatly improved compared to one-single BBO crystal. Figure 7(a) shows a single- 
stage compression simulation with a low soliton order. Figure 7(b) shows that high-order soliton 
compression may reach the limit of pulse duration (6 fs) but leave most of energy in the pedestal 
(Q^ = 0.23); while keeping low-order soliton compression in two-section BBO crystals can 
hold high-quality pulse (Qc = 0.72) with maximum peak intensity while minimum pulse width 
approaching the same duration as shown in Fig. 7(c). Figure 7(d) plots the pulse intensities at 
the optimal positions in one- and two-section BBO crystals. 

4. Conclusion 

We proposed a new scheme for compressing longer >> 100 fs pulses through a soliton effect 
to few-cycle duration with high quality. The cascaded quadratic nonlinearity was engineered 
through varying the phase mismatch in ways of multi-section quasi-phase-matching (QPM) 
structure with different domain periods and bonded BBO crystals with different orientation 
of the optical axis. By adjusting the residual or QPM phase mismatch, we made A^eff > 1 in 
the beginning of each section and kept soliton compression in the stationary regimes, so that 
clean few-cycle pulse compression was accessible. At the same time, each section was designed 
to keep a low effective soliton number, as soliton compression of pulses with small soliton 
numbers can achieve a small pedestal and high pulse quality. 

We in detail presented the numerical results of soliton compression in multi-section QPM 
structure at the wavelength of 1.56 /im. The available cascaded nonlinearities are compared in 
unpoled lithium niobate (LN) and periodically poled lithium niobate (PPLN), and we showed 
by using QPM the advantage of enough QPM phase-mismatch space for second- stage compres- 
sion. However, the shortest available wavelength for pulse compression in the QPM scheme is 
around 1.37 /im, due to competing Kerr nonlinearities. The wavelength cut off is substantially 
higher than that in LN, where compression may occur even down to around 1 micron. Examples 
of soliton compression to 13 fs (below three optical cycles) were summarized in multi-section 
QPM with different input pulse duration and soliton numbers. Input pulses with large soliton 
number experienced a stronger competing Kerr XPM effect, which is self-focusing and there- 
fore reduced the overall self-defocusing nonlinearity. 

Numerical results showed that the compressed pulse with less than three-cycle duration was 
close to the limit of the higher-order soliton compression while keeping quality factor more 
than 0.40. We also discussed that the tolerance of input intensities and the error of the QPM 
structures affect the soliton compression in the multi-section QPM. The Raman contribution to 
few-cycle soliton compression was also discussed. 

Finally we implemented the scheme in bonded BBO crystals at the wavelength of 1.03 /im. 
The numerical results confirmed that pulse quality of soliton compression in two-section BBO 
crystals is greatly improved while minimum pulse width approaching the same duration (6 fs) 
compared to one-single BBO crystal. 
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